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Abstract 

We prove locality estimates, in the form of Lieb-Robinson bounds, for classical oscillator 
systems defined on a lattice. Our results hold for the harmonic system and a variety of anhar- 
monic perturbations with long range interactions. The anharmonic estimates are applicable to 
a special class of observables, the Weyl functions, and the bounds which follow are not only 
independent of the volume but also the initial condition. 

1 Introduction 

A notion of locality is crucial in rigorously analyzing most physical systems. Typically, sets of 
local observables are associated with bounded regions of space, and one is interested in how 
these observables evolve dynamically with respect to the interactions governing the system. 
In relativistic theories, the evolution of a local observable remains local, i.e. the associated 
dynamics is restricted to a light cone. For non-rclativistic models, such as those we will be 
considering in the present work, the dynamics does not preserve locality in the sense that, 
generically, an observable initially chosen localized to a particular site is immediately evolved 
into an observable dependent on all sites of the system. 

In 1972, Lieb and Robinson [5] explored a quasi-locality of the dynamics corresponding to 
non-relativistic quantum spin systems. Roughly speaking, a quantum spin system is described 
by a self-adjoint Hamiltonian, which describes the inteactions of the system, and its associated 
Heisenberg dynamics, see e.g. [2] for more details. The estimates they proved, which we will 
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refer to as Lieb-Robinson bounds, demonstrate that, up to exponentially small errors, the thiie 
evolution of a local observable remains essentially supported in a ball of radius proportional to 
v\t\ for some v > 0. This quantity v defines a natural velocity of propagation, and it can be 
estimated in terms of the system's free parameters, for example, the interaction strength of the 
Hamiltonian. 

The models analyzed in this paper will correspond to a classical system of oscillators evolving 
according to a Hamiltonian dynamics. Hamiltonians of this type have frequently appeared in 
the literature as their analysis provides an important means of studying the emergence of non- 
equilibrium phenomena in macroscopic systems. For example, rigorous results on the existence 
of the thermodynamic limit for these models date back to [9]. A notion of quasi- locality, sim- 
ilar to the Lieb-Robinson bounds mentioned above, for these oscillator systems was originally 
considered in 1978 by Marchioro et. al. in 15j, and a recent generalization of these estimates 
appeared in [4] . Both of these results were obtained in the spirit of deriving an analogue of the 
Lieb-Robinson bounds found in f8]. 

Over the past few years a number of important improvements on the original Lieb-Robinson 
bounds have appeared in the literature [12l ISl [lOl |TT] , see [M] for the most current review 
article. These new estimates have found a variety of intriguing applications [SI [31 [T31 [3 1 but 
perhaps most interestingly for the present work, the results found in |11| establish bounds 
which are applicable beyond the context of quantum spin systems. In [11) . the authors prove 
a version of the Lieb-Robinson bounds for quantum anharmonic lattice systems. Motivated by 
the methodology introduced in [TT] , we are returning to the classical setting to re-derive distinct 
estimates for anharmonic lattice systems. 

To express our locality results more precisely, we introduce the following notation. We will 
consider systems confined to a large but finite subset A C Z'^; here i/ > 1 is an integer. With 
each site x G A, we will associate an oscillator with coordinate £ R and momentum £ R. 
The state of the system in A will be described by a sequence x = {{Qx,Px)}xeA, and phase space, 
i.e. the set of all such sequences, will be denoted by Xa- 

A Hamiltonian, H, is a real-valued function on phase space. Typically the Hamiltonian of 
interest generates a fiow, $t, on phase space. Specifically, given H : Xa R one defines, for 
any i e R, a function $t : — > by setting $t(x) = {{qx{t),Px{t))}x£A, the sequence whose 
components satisfy Hamilton's equations: for any x G A, 

. , , dH , , , , 

qx{t) = ^($t(x)), 

. , . dH , ^ , , , ^ 
oqx 

with initial condition $o(x) = x. 

To measure the effects of this Hamiltonian dynamics on the system, one introduces observ- 
ables. An observable A is a complex-valued function of phase space. We will denote by Aa the 
space of all local observables in A, i.e. the set of all functions A : X\ — > C. A given Hamiltonian, 
H, generates a dynamics at on the space of local observables in the sense that, for any i e R, 
the dynamics at : Aa is defined by setting at{A) = Ao^f 

For the locality result we will present, the notion of support of a local observable is important. 
Given A G Aa, the support of A is defined to be the minimal set A C A for which A depends 
only on those parameters Qx or with x £ X. 

As in [15] , see also |4] , our locality result will be expressed in terms of the Poisson bracket 
between local observables. Here the Poisson bracket is the observable given by 

{AB}-y"— — - — — (12) 
leA '^P^ ' 

for sufficiently smooth observables A and B. 
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Observe that for disjoint subsets X, F C A and observables A with support in X and B 
with support in Y, it is clear that {A, B} = 0. The quasi-locahty question of interest in this 
context is: given a Hamiltonian H, its corresponding dynamics at, and a pair of observables 
A and B with disjoint supports, is there a bound on the quantity {at{A),B} for small times 
t7 Physically, one expects that if the Hamiltonian is comprised of local interaction terms, then 
there should be a bound on the velocity of propagation through the system. Such intuition 
could be confirmed by establishing an estimate of the form 



I {at (A), 5} (x) I < Ce-^('^(^'^)-''l*l), 



(1.3) 



where d{X,Y) denotes the distance between the supports of the local observables A and B. 
This bound demonstrates that for times t with |t| < d{X,Y)/v, the Poisson bracket remains 
exponentially small, and the number v > appearing above is a bound on the system's velocity. 
In proving estimates of the form (|1.3p . special attention must be given to the dependence of 
the constants C, fi, and v on the observables A and B, the initial condition x, and the free 
parameters in the Hamiltonian. Most crucially, these constants must be independent of the 
underlying volume A, so that they persist in the thermodynamic limit; once the existence of 
such a limit has been established. 

As we have mentioned before, bounds of the form p.3p have appeared in the literature, see 
|15j and more recently [4], for a variety of different Hamiltonians. Both our approach and our 
estimates are distinct from those mentioned above. For example, we do not work with time 
invariant states, and our bounds are independent of the initial conditions. Our main goal is to 
provide a new method for establishing these estimates, and we are strongly motivated by the 
quantum techniques found in [llj . 

We begin by considering finite volume restrictions of the harmonic Hamiltonian, i.e. 
A:a — > K is given by 



xeA 



(1.4) 



where e^, for j = 1, . . . , j/, are the cannonical basis vectors in Z", and the parameters w > and 
Xj > are the on-site and coupling strength, respectively. As is well-known, a variety of explicit 
calculations may be performed for this harmonic Hamiltonian. Perhaps most importantly, for 



any integer L > 1 and each subset A^ = {—L,L]'^ C Z'', the flow ^^'^^ corresponding to Hj^ 



may be explicitly computed, see Section 12. II for details. Once this is known, a locality estimate 
easily follows for a specific set of observables. We will equip the set of local observables Aa^ 
with the sup-norm, and we will say that A G Aal is bounded if 



\A\ 



sup 



l^(x)| 



(1.5) 



(1) 



is finite. Furthermore, we will denote by 
xeAL, §^ e AAr., §^ e Aa^,, and 



the set of all A G for which: given any 



sup max 

xGAl 





dA 




dA 




( 


dqx 


1 

oo 


dpx 


OO-^ 



< oo . 



(1.6) 



We can now state our first result. 



Theorem 1. Let X and Y he finite subsets of 1/ and take Lq to be the smallest integer such 

For 



that X,Y C Aif, . For any L > La, denote by a^'^ the dynamics corresponding to Lf^'- 



any fi > and any observables A, B £ -^Al support of A in X and support of B in Y , 

there exist positive numbers C and vt, both independent of L, such that the bound 



|a^^(A),i3| <C\\dA\\^\\dB\\^mm{\X\,\Y\ 



,-l-i(d{XX)-Vi,\t\) 



(1.7) 
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holds for all t CzM. 

Some additional comments are in order. First, the quantity d{X, Y) appearing above denotes 
the distance between the sets X and Y, measured in the L^-sense, and for any Z C A^, the 
number \Z\ is the cardinahty of Z. Next, the fact that the bound (|1.7p is true for any /i > 
imphes that the Poisson bracket above has arbitrarily fast exponential decay in space. To 
achieve faster decay in space, however, the numbers C and Vh increase. We describe an optimal 
harmonic velocity Vh{fJ-) in Section [2.21 

One novelty of our approach is that the bound in (I1.7P is not only independent of the length 
scale L, it is also independent of the initial condition x £ Xj^^. This fact remains true when we 
consider anharmonic perturbations, see Theorem [2] below, and thereby distinguishes our result 
from that of [TS] and [1]. 

Our next result, Theorem[2]below, concerns single site perturbations of the harmonic Hamil- 
tonian. To state this precisely, fix a function y : K ^ M. For any site z £7L^ define Vz : Xt^,^ — > R 
by setting Vz{^ = l^(9z)- We consider finite volume anharmonic Hamiltonians H'^^ : Xaj^ M 
of the form 

In order to prove Theorem^ we need the following assumptions onV: V G C^(R), V G L^(R), 
V" e L°°(R), and 



V'{r) 



dr < CX3 , 



Here V is the Fourier transform of V . Under these assumptions, we prove a locality result 
analogous to Theorem [T] As is discussed in Section [231 see also the proof in Section [31 a specific 
class of observables, the Weyl functions, are particularly well-suited for our considerations, and 
they are defined as follows. For any function / : A^ C, the Weyl function generated by /, 
denoted by W{f), is the observable W{f) : Xa^ C given by 



[VF(/)](x) = exp 



Re[fix)]q,+lm[fix)]p, 



xeAi 



(1.9) 



Clearly, if / is supported in AT C A^, then W{f) is supported in X as well. Moreover, it is easy 
to see that for any function / : Al C, ||M^(/)||oc- = 1- Our next result is 

Theorem 2. Let V -.R-^R satisfy V € C^iR), V € L^(R), V" e L°^{R), and ny , as defined 
above, is finite. Take X and Y to he finite subsets of and let Lq to be the smallest integer 
such that X,Y C A^q. For any L> Lq, denote by the dynamics corresponding to H^'^ . For 
any /x > and any functions /, g : A^^ — > C with support of f in X and support of g in Y , there 
exist positive numbers C and Vah, both independent of L, such that the bound 

||K(T4^(/)),^y(5)}L<C||/|U||5lloomin(|X|,|y|)e-^('^(^^^)-''- (1.10) 

holds for all t CzR. 

The assumptions on V above are sufficient to imply that V is bounded. For this reason, our 
results do not apply to more substantial perturbations, e.g., those of the form Vz{x) = q^. In [4], 
the authors do consider, for example, quartic perturbations. They prove that, for reasonable time 
invariant states and a set of initial conditions of full measure, after a time t local perturbations 
of thermal equilibrium are exponentially small in log^ (t) at a distance larger that t log" (t) . This 
is insufficient to conclude the existence of a finite velocity u > as we have discussed above. It 
is an interesting question to determine whether or not this genuinely describes the behavior of 
such systems. We will not answer this question in the present work. 
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The paper is organized as follows. In Section [21 we discuss our results concerning the Har- 
monic Hamiltonian and prove Theoreni[TJ Using an interpolation argument, we prove Theorem[2] 
in Section [3l This result demonstrates that our locality bounds for Weyl functions are preserved 
under certain single-site anharmonic perturbations. In Section lU we generalize Theorem [2] to 
cover a wide class of multi-site perturbations. Finally, Section [5] contains a variety of useful 
solution estimates used throughout the paper. 

2 The Harmonic Hamiltonian 

The main goal of this section is to prove Theorem[TJ For the convenience of the reader, we begin 
with a subsection describing some basic features of the harmonic Hamiltonian. In particular, 
we reintroduce the Hamiltonian and find an explicit expression for the corresponding flow. In 
the subsections which follow, we prove two locality estimates. The first is valid for a general 
class of smooth and bounded observables. The next holds for a special class of observables, the 
Weyl functions. This latter result will be particularly useful in subsequent sections. 

2.1 Some basics 

For any integer i > 1, we consider subsets Kl — {~L, LY C 2^ and the finite volume harmonic 
Hamiltonian -ff;^^ : Xa^ K given by 

i^^(x)= ^ p2^c.2g2 + ^A,(g,-fc+e,)'. (2.1) 

aJSAi, j = \ 

Here, for each j = 1, . . . , z^, the are the canonical basis vectors in Z", w > 0, and \j > 0. The 
model in (j2.ip is defined with periodic boundary conditions, in the sense that Qx+cj ~ Qx-i2L-i)ej 
if a; G Ai but x + Cj ^ . 

Our first task is to provide an explicit expression for the flow corresponding to (|2.ip . In 
doing so, we will fix an integer value of i > 1 and drop its dependence in a variety of quantities 
to ease notation. Given any x S X^j^ and t e M, the components of $J'(x) = {{qx{t),Px{t))}xeAL 
satisfy the following coupled system of differential equations: for each x G Al and t G M, 

qx{t) = 2px{t) 

Px{t) - -2Lu\x{t)-2j2^j{'^'lAt)-1x+e,{t)-qx-e,{t)) (2.2) 

i=i 

with initial condition {{qx{0),Px{0))}xeAL — ^- Introducing Fourier variables, the system de- 
fined by (|2.2p decouples which leads to an exact solution. This is the content of Lemma [3] found 
below. 

Before stating Lemma [3l it is useful to introduce some additional notation. Fourier sums 
will be defined via the set A^ given by 

Note that A^ C (-tt, tt]'^ and |A^| = \Al\ = {2LY . The following functions play an important 
role in our calculations. Suppose lo > Q and take 7 : A^ ^ K. to be given by 



w2 + 4^AjSin2(fcj/2), (2.3) 
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and for each m G { — 1, 0, 1} and any t G R, set h^' : ^ M to be 



h[ (.x) ~ Im 
hf\x) = Re 



1 ^ gi(fc-2;-27(fe)t) 



E 



— y 



(2.4) 



|A 



feeA* 



Each of these functions depend on the length scale L, however, we are suppressing that depen- 
dence. 

Lemma 3. Suppose w > 0. For any x G X/^^ and t G K, the mapping $^ : Aaj^ — + Aa^^ 
is well- defined. In particular, for each x G and t G M, </ie components of <&J'(x) = 
{(9a;(t),Pa;(i))}^gAi, "'"^ given by 



(2.5) 



(2.6) 



E %{^)h^t\^-y)+Py{^)h^t\^-y)- 

Here, if necessary, the function values h["^^ {x — y) are defined by periodic extension, and we 
regard x= {{qx(0),Px{0))}xeAL ■ 

Proof. Taking a second derivative of (|2.2p . we find that for each x G A^ and any t G M, 



(2.7) 



For any fc G A^ and i G M, set 



(t) - E and Pfe(^) = '^''"'"P^t)- (2 



xeAi 



Inserting p.7p into the second derivative of (|2.8p . we find an equivalent system of uncoupled 
differential equations. In fact, for each k G A^ and any i G K, 



Qk{t) = -4io^Qk{t) - 4^ A, (2 - e^'^^ - e-^^^) Qu{t) = ^^i{kfQk{t) 
Pkit) = -4c^2p^_(^) _ 4^ A, (2 - e*^-^ - e"*^'^) Pfc(t) = ~Aj{k)^Pk{t) , 

where 7 is as in (|2.3p . The solution of (|2.9p is given by 

Qfc(t) =Cfee-2^^W*+^e2*^('=)* 



(2.9) 



Pfc(i) = Dke 



D-k.e 



2i^(k)t 



(2.10) 
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where — fc is defined to be the element of whose components are given by 

^ 1 ""i otherwise. 

The relationship between the coefficients in (|2.10|) above is derived using the fact that the initial 
condition is real- valued, e.g., 



Qfc(0)-Q-fe(0) and gfc(0) = g_fc(0). 
Using Fourier inversion, we recover the components of the ffow from (|2.10p . In fact, 

qS) - E 

= y C7fce*(*'-"-2''('=)*' +Cfce-^('=-"-2''('='*), (2.11) 

and similarly, we find that 

p^it) = jj^^r(k-x-2-^(k)t) ^^^-r(k-x-2-^(k)t) _ (2.12) 

^\k^l 

To express these solutions explicitly in terms of the initial condition, we observe that 



Qfc(0) = Cfe + C_fc and Pk{Q)=Dt, + D^k, (2.13) 

and introduce 



Bk = -jL=pm~^f-^Qk{^) with . -jL=p_,m^f-^Q-k(^)- (2.14) 

It is easy to see that 



Qk{Q) = -1== {Bk - B^k) and Pfe(O) = \l ^ {Bk + B^k) , (2.15) 
\/27(fe) 



and therefore, 



C,^^£^ and Dk^xh^Bk. (2.16) 
^2^ ^ Y 2 ^ ' ' 

Plugging this into (|2.1ip . we find that 

qJt) ^ 1 ^g^ ^»(fc..-27(/c)t) _ iBk ^_,^k-x-2^ik)t) (2.17) 

VWTlk'klV^) N/27(fc) 

= I y gfc(0)e'('=-"-27(fe)t)+g^g-»(fc-x-27(fc)t) 



-Pfc(O) ^j(fc.:r-27(fc)t) _ Pk{0) ^-i(k-x-2i(k)t) 



E 



7^ E [g.(o)e^(--^('=)*)] ^ 7^ E I- 



-Pfc(O) ^i(k-x-2j(k)t) 



7(fc) 
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Moreover, one finds that 



Re 



Ofc(0)e 



i{k-x~2'y{k)t) 



yeAL 



while 



Im 



Pk{0) i(k-x-2',(k)t) 



7W S 



^t{k-ix-y)-2^ik)t) 



J(k-(x-v)-2-,{k)t) 



7(fc) 



: Qo riofinori (!'; /Ill conclude that 



With the functions hi , as defined in 



(2.18) 



(2.19) 



(2.20) 



as claimed in (|2.5p . A similar calculation yields (|2.6p . Since the functions are real valued, 
so too are the solutions qx (t) and Px {t) . This proves Lemma [31 □ 



Remark 4. analogue of \2. 5]) and h2. b]) holds in the event that lo — 0. This is seen by 
proceeding as in the proof of Lemma\^ and observing that now 7(0) = 0, but "f{k) ^ for k ^ 0. 
For k ^ 0, the formulas above are correct, and a simple calculation shows that, in this case, 



Qoit) = Qo(0) + 2Po(0)t 
Po{t) = Po(0), 



(2.21) 



similar to i2.10\) . One easily sees that the equations h2. 5\) and V2. b]) still hold with the convention 
that 

J{k-x~2y{k)t) 



h[-'\x) 



2t 



Im 



1 



|Al| 



E 



7(fc) 



(2.22) 



fc6A*\{0} 

We end this subsection with the following crucial estimate which was proven in [11 . 

Lemma 5. Fix L > 1 and consider the functions /ij™'' as defined in |^.^[ ) for m G { — 1,0,1}. 
For any /i > 0, the bounds 

^ g-M(b|-c.,.max(^,e<^'/^) + i)|t|) 



,(-1) 



(x) 



h^^\x) 



< ^-l^g-p(k|-c..^max(f ,e(^/^) + i)|t[) 

< c . e^/2e-^(l-l-^-^ • -'^^''+')l*l) 



(2.23) 



hold for all t M. and x £ A^. Here \x\ = X^J^i ^'^^ ""^^ may take c^^ \ = (iJ^ + 

We refer the interested reader to Lemma 3.7 of [11' for the proof. Moreover, we stress that 
Lemma [S] is valid for all a; > 0. 



2.2 A general locality estimate 

Our first locality bound for the harmonic Hamiltonian follows directly from Lemma [3] and 
Lemma O We state this as Theorem [HI below. As we will see. Theorem [T] is an immediate 
consequence of Theorem |6l Recall that we have defined J^]^^^ to be the set of observables 
A G Ak^ for which: given any a; G Al, S -Aa^, ^ G Ak^, and 



\\dA\\ 



sup max 





dA 




dA 




( 


dqx 


1 

OO 


dpx 





< OO . 



(2.24) 
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Theorem 6. Let X and Y be finite subsets of 'L^ and take Lq to be the smallest integer such 



that X,Y CZ A^g. For any L > Lq, let a^ '^ denote the dynamics corresponding to Hj^'' . For 
any /i > and any observables A, B E -^aI with support of A in X and support of B in Y , the 
bound 



|af'^(A),s} <C\\dA\\^\\dB\l^ J2 e~'^(''(^'«)-"-^'°^''(M'^*'''''^')l*l) (2.25) 
holds for all t G R. Here 



xex.yeY 



dix, y)=y] min \xj - y, + 2L r]j \ 



(2.26) 



is the distance on the torus and the constants may be taken as C = {2 + c^j.ac^/^ + c^\) with 
c.,a-(c.^ + 4E;.iA,)V2. 



Proof. The Poisson bracket is easy to calculate. In fact, for any x G , 



{of-(A,.«}](x,= i:AA(*nx, 



dB 



d 



dpy ^ dpy 



(*f'(x,).|£(x,. (2.27) 



dq. 



By the chain rule, 



OA 
dqx 



dqx 
dqy 



dA 
dpx 



dpx 



(2.28) 



and a similar formula holds for -^A ' (x) j . Now estimating (|2.27p . we find that 
{a,^'^(A),B}||^ < \\dA\U\dB\\^ \hi-'\x-y)\+2\hi'\x-y)\ + \h['\x-y) 



xex,yeY 

using Lemma [3l The bound in p.25p now follows from Lemma [5l 



(2.29) 
□ 



From Theorem [6l and specifically the bound ()2.25p . we see that for any /i > 0, the harmonic 
velocity Vh is essentially described by 



«h(/^) = c^,Amaxf-,e(''/2)+i' 

In fact, given (|2.25p for some /i > 0, it is easy to see that for any < e < 1, 

J2 e-''^(^'«) < e-^'"*(^'^) min (|X|, ^ e-'^^^-^^'^^"'^) 



(2.30) 



(2.31) 



xGX,yGY 



where we have set d{X, Y) = miiixi=x,yeY d{x, y). Thus, Theorem [T] is a simple consequence of 
Theorem ini It is interesting to note that for any L the quantity 



(2.32) 



where \z\ denotes the L^-mctric on . Given this and the fact that, for sufhciently large L, 
the distance d{X, Y) agrees with the L^-distance between X and Y , it is clear that the estimate 
proven in Theorem [T] is genuinely independent of the length scale L. 
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Since the bounds are valid for any fi > 0, Thcorcm|n]demonstrates arbitrarily fast exponential 
decay in space with a velocity that depends on fi. Typically, however, one is interested in the 
best possible estimates on Vh given some decay rate. In this sense, the optimal harmonic velocity, 
as described by ()2.30|) . occurs when the equation 

I = e(''/2)+i (2.33) 

holds. It is easy to see that the solution to (|2.33|) . denoted by /io, satisfies 1/2 < /io < 1, and 
therefore the corresponding velocity Uh(/^o) < 4:Cuj,x- 

2.3 The harmonic evolution of Weyl functions 

In preparation for our arguments in Sections [3] and [H we will now present a different proof of our 
locality result, analogous to Theorem [SI valid for Weyl functions. Recall that a Weyl function 
is an observable, generated by a function / : — > C, with the form 



[iy(/)](x) = exp 



J2 Re[/(a;)]fe+Im[/(x)]p, 



(2.34) 



One important property of the Weyl functions is typically referred to as the Weyl relation. 
We state this as Proposition [71 

Proposition 7 (Weyl Relation). Let f,g : Al C. We have that 

{W{f), W{g)} = -Im[(/, g)]W{f)W{g). (2.35) 
where the inner produet is taken in 1'^{Kl)- 
Proof. A direct calculation yields 

{WU).W{g)} = ^W{f)^W{g)--^W{f)-^W{g) 

= (-Re[/(x)]Im[g(x)] + lu,[S{x)]Re[g{x)])W{f)W{g). 

Noting that 

Im[(/,5)] = Im[ ^ W)g{x)\ 

= (-Ini[/(x)]Re[.9(a;)] + Re[f(x)]lui[g{x)]) 

x&Al 

proves the proposition. □ 

Another useful property of the Weyl functions is that the harmonic dynamics leaves this 
class of observables invariant. This important fact, which follows immediately from Lemma [31 
is the content of the next proposition. Before stating this, it is convenient to introduce notation 
for the convolution of two functions /, (7 : Al ^ C, 

U*9){x) = Y fiy)9i^-y)^ (2-36) 

where, if necessary, g{x — y) is calculated by periodic extension. 
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Proposition 8. Let f : Al C and take t G R. 

a'^'^iWU))^Wift), (2.37) 
where 



/* = / * (^/^r + 2^^*^^^ + hn) + f * (/^^^^ - hr') i • (2.38) 

with h[ ^\h[^\ and h^^^^ as in ^2.4^ . 
Proof. For any point x G Aai, , we have that 

{W{f))]{x) = expjz ^ Re[/(x)]q,(t)+Im[/(x)]p,(t)] 



h,L/ 

a. 



exp|z ^ Re[f{x)] ^ qy{0) {x - y) - py{0) h[-'\x - y) 
x&Al y&^L 

xexp ji ^ Im[/(x)] qy{0)hi'\x~y)+py{0)h['\x-y) 
\ xeAL yeAL 

= exp [ z ^ qy{0) J2 Rc[f{x)]hi"\x -y)+ Im[/(x)] h['\x - y) 

\ yeAL xeAL 

xexp jz ^ Py{Q) ^ lm[f{x)]h[°\x-y)-Re[f{x)]h[''\x-y) 

\ yeAL xeAL 

= [W^(/t)](x), (2.39) 

where we have defined the function ft : — > C by (|2.38p . □ 

It is obvious that Theorem [9] below follows immediately from Theorem \6\ since the Weyl 
functions are clearly in A^j^j^ . We will here give a different, but equally short, proof which uses 
the specific properties of Weyl functions. 

Theorem 9. Let X and Y be finite subsets of and take Lq to be the smallest integer such 
that X,Y C ^Lo- For any fi > 0, L > Lq, and any functions f,g : A^^ — > C with support of f 
in X and support of g in Y , the hound 

|a^^(W^(/)),ly(5)| <C||/||oo||.g||oo V e-K'i(-.J/)------(^^™^)l*l) (2.40) 

xex,yeY 

holds for all t G M. Here, as in 112. 26\) . d{x, y) is the distance on the torus and the constants 
may be taken as C ^ {1 + c^^xe^^/"^ + cj\) with c^^\ = {lo'^ + 4X]j=i ^j Y^"^- 

Proof. Combining Propositions [8] and [71 it is clear that 



[a1-'\W{f)),W{g)] = {W{ft).W{g)} = -lm[{ft, g)]W{ft)W{g). (2.41) 
In this case, the bound 

|a^^M^(/)),Iy(ff)} <|Im[(/„5)]|, (2.42) 
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readily follows. Appealing again to Proposition [SJ we have that for any y £ 



Y^f(a:)('-hi'\x-y)-'-hi-'\x-y) 



(2.43) 



and therefore, 



\{fu9)\ = 



ft{y)9{y) 



< ||/||oo||g||co E \h["\x-y)\ + \h[-'\x~y)\ + \h['\x-y)\. (2.44) 

xeX.yeY 



Theorem |9] now follows from Lemma [3 



□ 



We end this section with a corollary of Theorem [9] that will be particularly useful in the next 
sections. The locality bound we prove for the anharmonic dynamics is derived by iterating a 
certain inequality involving the harmonic estimate. With this in mind, it is useful to introduce 
the following family of decaying functions. For any ^ > 0, consider : [0, oo) — » (0, oo) defined 

by 

(2.45) 

Clearly, these function F^^ also depend on the quantity v > 1, which is the dimension of the 
underlying lattice in our models, but we will suppress that dependence in our notation. Unlike 
the bare exponential e~^^ , these functions have the following nice property. There exists a 
number Ci/ > for which, given any pair of sites x, y G Z"^, 



^ F,{\x - z\)F^{\z - y\) < CM\^ - y\) ■ 



(2.46) 



Here one may take 



C, 



= 2"+! V 1 

(l + |z|)''+i 



(2.47) 



Functions of this type were introduced in [10] , see also [llj , as an aide in proving Lieb- Robinson 
bounds. We will use them here as well. 

We can rewrite the decay expressed in our harmonic estimates, i.e. (|2.25[) . in terms of these 
functions F^. 

Corollary 10. Let X and Y be finite subsets of II' and take Lq to be the smallest integer such 
that X,Y C A^,-, . For any fi > 0, e > 0, L > Lq, and any functions /, g : A^^ — )■ C with support 
of f in X and support of g in Y , the bound 



{a,'''^(iy(/)),M/(.9)}L<C||/|U||5||ooe('^+^)"'-(^+^)l*l FMx,y)) , (2.48) 

holds for all t (zR. Here 



C 



(1 + c^,xe^ + cZ\) sup [(1 + sY+^e-^^] 

s>0 



(2.49) 



and v\i is as defined in i2.30\) . 
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3 Single Site Anharmonicities 



In this section, we will prove a locality result, analogous to Theorem [SI for a specific class of 
perturbations of the harmonic Hamiltonian. A much more general result, which follows from 
the same basic arguments, is presented in the next section. We begin with a precise statement 
of the models we consider, and then prove the result. 

To make our basic technique more transparent, we will only consider single-site potentials 
that are generated by a particular function V in this section, see Section S] for a more general 
result. Let y : E ^ M satisfy V G C^{R), V G V" G L°°(R), and suppose further that 

the quantity 

,1V = f \r\ \V'ir)\dr (3.1) 
Jr 

is finite. Here V is the Fourier transform of V' . Given such a function V and an integer L > 1, 
we define an anharmonic Hamiltonian H^^ : Aa^, ^ K by setting 

jjAl =Ht-+Y,y^^ (3-2) 

z6Ai 

where for each z G A^, the function Vz : -^Ai, ^ R is given by K(x) = V{qz)- 

As is discussed at the end of Section 12.31 we will state our result in terms of the functions 
: [0, oo) (0, oo) given by 

with v > corresponding to the dimension of the underlying lattice . The goal of this section 
is to prove the following result. 

Theorem 11. Suppose V : R ^ R satisfies V G C^iR), V G i^(M), V" G L°°{R), and 
Kv , as in 113. 1]) above, is finite. Let X and Y be finite subsets of and take Lq to be the 
smallest integer such that X,Y C A^^. For any L > Lq and t G M, let denote the dynamics 
corresponding to H^'^ . Then, for any /i > 0, e > 0, and any functions f,g : A^^ — > C with 
support of f in X and support of g in Y , the bound 

||{af(M/(/)),M^(g)}||^<C|l/|U||5llooe*l*' ^ F^(d(x,y)) (3.4) 

x£X,y£Y 

holds for all t £R. Here one may take 

C = (1 + c^,xe^ + cZ\) sup [(1 + sr+^e~''] (3.5) 

and 

6 = (5(/i, e) = (/i + e)vi,{fi + e) + CC^/ty (3.6) 

where Vh is as in 112. 30\} . Cu is in |^.^7| j, and kv is in 113.1]) . 

Before we prove Theorem lll[ we comment on the corresponding anharmonic velocity. With 
arguments similar to those given after the proof of Theorem [6l it is clear that Theorem [11] 
implies Theorem [2] In this case, we find that an upper bound on the anharmonic velocity for 
this model is 

i'ah(M,e)= fl + -Vh(A' + ^) + ^^^- (3-7) 

We now present the proof. 
Proof of Theorem ] 111 Our proof of this estimate is perturbative, and we begin by interpolating 
between the harmonic and anharmonic dynamics. Fix L > Lq as in the statement of the theorem. 
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Since we will regard both the harmonic and anharmonic dynamics on the same volume A/,, we 
drop the dependence of each on L. Observe that for any t > 0, 

{atiWif)), W{g)} - {a'liWif)), W{g)} - ^ {«. K-.(W^(/))) , Wig)} ds . (3.8) 



Moreover, a direct calculation shows that 



(at, (Wif))) = a. ({at,(W^(/)), H}) - a, (at, {{Wif), H„})) 
= as{{al^{Wif)),H-H,,}) 

= ^ a.({at,(W^(/)),V;}) . (3.9) 

The Poisson bracket on the right-hand side of (|3.9p can be simplified 

{ats{W{f)),V,} = {W{ft-s),V,} = ~tlm[ft-s{z)] W{ft-s)V^. (3.10) 
For the first equality above we used Proposition [51 and we have denoted by the function 

-.Xa.^r with v;'(x) = v'iq,). 

These calculations lead to a particularly simple differential equation and thus, eventually, 
the bound p.l7p appearing below. In fact, for t > fixed and < s < t, define the function 

^t{s) = {aM-siW{m,W{g)}. (3.11) 

We have shown that 

^Ms) = ^ {a.({at,(M/(/)),K}),m5)} 

zeAt 

- i/:t(s)^t(s) + Qt(s), (3.12) 



where 



Ctis)^- J2 Ini[/t-.(z)]a.(K'), 

Qtis) = E K-.(W^(/))) {a.(K'),W^(5)} 



(3.13) 



and the final equality in p.l2|) follows from an application of the Leibnitz rule for Poisson 
brackets. Since for each fixed s, Ct{s) is a real- valued function of phase space, the solution Ut 
of 

^Utis) = ~iCtis)Utis) with Ut{0) = 1, (3.14) 
ds 

is a complex exponential. In addition, it is easy to see that 

£(^^,,{s)Ut{.s)) ^ Qt{s)Ut{.s), (3.15) 

and therefore. 



^iJt{t)Ut{t) = ^t{0) + / Qtis)Ut{s)ds, (3.16) 
Jo 

from which 

\\{atiWif)),Wigm^ < ||K(W^(/)),W^(.9)}|L 

J2 f |Im[/t-.(^)]| \\{as{V^),W{g)}\\^ ds , (3.17) 

,= A. "'0 



zeAi 
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readily follows. 

Now, if Vj' was a Weyl function, then we could immediately iterate the inequality in (|3.17[) 
and derive a bound. This is not the case, however, our assumptions on V allow us to write 
as an average of Weyl functions through its Fourier representation. In fact, we write the Fourier 
transform of V' as 

V'ir)^^ f e-^'^^V'{q)dq, (3.18) 

and by inversion, one has that 

V'iq) = [ e'''«T7'(r)dr. (3.19) 



This implies that, as a function of phase space, can be expressed as 

T/' = / W{rS^)V'{r) dr (3.20) 



where rSz '■ ^ M is the function that has value r at z and otherwise. Inserting p.20p into 
(|3.17p . we have that 

\\{atiWif)),Wigm^ < \\{a1iWif)),W{g)}\\^ (3.21) 



J2 f \l^[ft-s{z)]\ I 



V'{r) II {as {WirS,)) , W{g)} |L dr ds . 



At this stage, we can finally iterate the inequality. First, however, we insert the harmonic bound 
found in Corollary (TUl 

Recall that for any /i > and e > we have established (|2.48p with a constant C as in (|3.5p . 
With equation (|2.43p . it is easy to see that, for any /i > and e > 

\lm[Mz)]\ < Cl|/l|ooe(^+^)'"'(^+^)l*l J2 F,id{x,^)), (3.22) 

also holds for any z G A^ and t G R. To ease the notation a bit, we will denote by w = 
(/i + e)t;h(/i + e). Using these bounds, the inequality in (j3.2ip now takes the form 

\\{at{W{f)), W{g)]\\^ < C ll/IU biloo e'' J2 ^'^(^(^' 

xeX.yGY 

l-t . (3.23) 

+ C\\f\\ooY. E Wx,z)) / / \V'{r)\\\{as{W{rS.)),W{g)}\\^drds. 

Upon iterating (|3.23p m > 1 times, we find that 

m 

\\{at{W{f)),W{g)}\\^<C\\f\\^\\gUe'' ^ ^ a„(x, t) + i?™+i(0, (3-24) 

x£X,y£Y n=0 

where 

aoix,y;t)=F^,id{x,y)), (3.25) 

ai{x,y;t) = Ct j \r\ \v'{r)\ dr ^ F^{d{x,z)) F^{d{z,y)) 

< CKvC,tF^id{x,y)) , (3.26) 
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and in general, 



an{x,y;t) 



n 



(fi f \rk\\v^{rk)\drk] ^ {d{x, z^)) ■ ■ ■ {d{zn,y)) 



< ^^^^^^^F,idix,y)), (3.27) 



for any 1 < n < m. In (j3.26p and ()3.27|) . we have used (|2.46p several times. 

From Lemma [171 found in Section O it is easy to see that the apriori estimate 

\\{as(W{rS,)),Wig)}\\^ < Ci |y | |r| ||g|U exp (C2 t^) (3.28) 

holds for all < s < t. Thus, for t > fixed, the remainder term Rm+i{t) converges to zero as 
m — > 00. In fact, 

i?™+i(t) < C, \X\ \Y\ ||/||oo||g|Ue^*+^-*^ ^^7^f"^\ (3.29) 

(to + 1)! 

We have proven that 

||{a,(W^(/)),M/(5)}|L <C||/|U||5||ooe('^+^'=-^'^)* ^ F,{d{^,y)) , (3.30) 
i.e. p.4p as claimed. □ 



4 Multiple Site Anharmonicities 

In this section, we will generalize Theorem[TT]in such a way that it covers perturbations involving 
long range interactions. As in the previous sections, we will be fixing some integer L > 1 and 
considering only finite volumes Kl C 1^ . 

We will introduce these perturbations quite generally and then discuss the assumptions 
necessary to prove our locality result. To each subset Z C A^, we will assign a function 
V{-; Z) : M. and a corresponding function of phase space Vz '■ R defined by setting 

Vzix) = V{{q,},ez;Z) . (4.1) 

Here {qz}z£Z is regarded as a vector in and the number V {{qz}zez', Z) is calculated by 
evaluating V{-;Z) with as the value in the z-th component for each z ^ Z. With this 
understanding, we will use the notation 

dzyz^-^Vz^dzV{-\Z), (4.2) 
oqz 

to denote the partial derivatives of Vz ■ 

In general, the finite volume anharmonic Hamiltonians we consider are of the form H^^ : 
R with 

^At = ^At + ^ 1/^ , (4.3) 

2CAt 

where the sum above is over all subsets of Kj^ . As we saw in Section [31 in order to prove our 
locality result, we need some assumptions on the functions Vz- We will now list these explicitly 
below. 

First, we use Lemma [TSl proven in Section [5l to provide explicit bounds on the components 
of the flow which, in particular, prevent the solutions from blowing-up in finite time. For these 
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estimates, we assume the perturbation above satisfies: 

i) For each Z C A^, the function Vz has well-defined first order partial derivatives. 

ii) There exist numbers C\ > 0, Ci > 0, and > such that for each x £ and any x £ 

2 

1 : ~ (4-4) 



zcAl / yeAL 



The decaying functions are as defined at the end of Section 12.31 

Next, much like in the proof of Theorem 111! we will need an apriori estimate on the Poisson 
bracket of specific, dynamically evolved observables. This is the content of Lemma [TTl found in 
the next section. To prove it we use Lemma [T51 and therefore, we must assume 

iii) For each Z C A^, the function Vz has well-defined second order partial derivatives. 

iv) There exist numbers C2 > and /i2 > for which: given any pair x,y G A^, the bound 

J2 \[d.dyVz] (x)| < C^F^, (d(x,y)) , (4.5) 

holds for all points x € . 

Lastly, we need the quantities that arise in our iteration scheme to be well-defined. For this 
we assume 

v) For each Z C A^, the first order partial derivatives of Vz are integrable. By this we mean 
that given Z C Al and z G Z, the function dzV{-;Z) is in L^(R^) with respect to Lebesgue 
measure. In this case, the Fourier transform of these functions exists, and we will write 

d^ir; Z) = —1^ / e— « d,V{q- Z) dq , (4.6) 



(27r)l^l 
for any r e R^. 

vi) For each Z C A^, we assume that the Fourier inversion formula holds for all first order 
partial derivatives of Vz- Thus, for any q G K^, 



d,V{q- Z) = J^^ e'^-" d,V{r- Z) dr , (4.7) 

and therefore, we will write 

[d.Vz] (x) = j^^ [W{t ■ 5z)] (x) d^{r- Z) dr , (4.8) 
where the function r ■ (5^ : A^ ^ R is given by 



[r ■ 5z]{x) = ^ 5z(x) hence \W{r ■ 5z)] (x) — exp 



(4.9) 



as required. 

vii) There exist numbers C3 > and /is > such that given any points x^y G A^, the bound 
V / \r\- VV{r-Z) dr <C^F^^{d{x,y)) . (4.10) 

x,yiaZ 

Here the vector-valued function W{-\ Z) : R^ has components dzV{-; Z) for each z G Z. 

The number |r|, corresponding to some r G R'^, is taken as \r\ = X^zez l^^l' seen in 

the proof below, any norm on R'^ satisfying \rz\ < ||r|| will suffice. As will become apparent 
below, we interpret the function F^^ in assumption vii) as our crucial estimate on the range of 
the interactions. 

We now state our most general result. 
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Theorem 12. Let X andY he finite subsets ojV and take Lq > 1 to be the smallest integer such 
that X,Y CZ A^f, . For any L > Lq and t € M., let denote the dynamics corresponding to the 
anharmonic Hamiltonian H^'^ in and suppose that the perturbation satisfies assumptions 

i) - vii) listed above. Then, for each e > and any functions f,g: Al,, C with the support of 
f in X and the support of g inY, 

\\{ai:{Wif)),Wig)}\\^<C\\f\U\g\\ooe'\'\ PM^^v)) (4-11) 

x£X.yeY 

holds for any t G K. Here 

C = (1 + c^.xe^ + c-\) sup [(1 + sY+'e-'^] . (4.12) 

and 

5^S{e) = {fi3 + e)vi,{p3 + e)+CC3Ci, (4.13) 

where Wh is as in Ii2.30\) . Ci, is in ^KJ^, and C3 is in ^.10\ ). 

One important difference between the bound we prove in Theorem 1121 above, in contrast to 
the one proven in Theorem [Til is that the spatial decay rate in (|4.1ip can be no greater than the 
rate /is appearing in (|4.10p . If fi^ > 0, then there is a corresponding velocity for this anharmonic 
system 

«ah(e)- fl + -Vh(M3 + e) + ^^^^^- (4.14) 

Since the case of /is = represents only polynomial decay in the interaction range, as measured 
by (|4.10p . the bound in (|4.1ip at most decays polynomially in distance between the supports of 
/ and g as well. 

Example 13. To clarify the general assumptions on the perturbation introduced above, we will 
consider a simple model with pair interactions generated by a single function. One can compare 
this example with the single site, anharmonic Hamiltonian analyzed in Section\^ Let 1/ : — » M 
be given and fix some number fi > 0. For each L > 1 and any Z d A^, define 

y(^,.^x(Pl^(.d{zi,Z2))V{-) ifZ={zi,Z2}, 

\ otherwise, 

and thereby, the anharmonic Hamiltonian 

H^^=Ht+ n^i,-.}' (4-16) 

zi,z2eAi, 

with V{zj_22}(x) — Ff^ {d{zi,Z2)) ■ V{qz^,qz2)- As one can easily check, the basic assumptions i) 
- iv) follow ifV has well-defined, second order partial derivatives and there exist numbers Ci, 
Ci, and C2 such that 

m&^\d,V{x,y)\<cA\x\ + \y\+Ci) (4.17) 

and 

max \d^djV{x,y)\<C2. (4.18) 
ije{i:2} 

// both first order partial derivatives of V are integrable and satisfy the Fourier inversion 
formula, then the condition vii) is satisfied when 

{\x\ + \y\) {\^d^{x,y)\ + \d2V{x,y)\^ dxAy < cx) . (4.19) 



Thus, under the above conditions, the model described by ( 
Theorem M^ and hence the corresponding locality result {4^.11 



^. J6p satisfies the assumptions of 



) is valid. 



18 



Proof of Theorem Much of the argument in the proof of Theorem [TT] also apphes here. 

Again, we fix L, regard both Hamiltonians on the same volume, drop the dependence of each of 
the dynamics on L, and interpolate. Let i > and set 

Ms) = {a.(at,(W^(/))), Wig)} (4.20) 

for < s < t. The calculation 

±asialAWif))) - 5] a,(K_,(iy(/)),^z}) 

ZcAl 

= ^ 5]lm[/,_,(z)]a,(W^(/,_,)).a,(9,Fz) , (4.21) 

ZcAl z£Z 

follows readily, and therefore, we derive a differential equation analogous to (I3.12p : namely 

-^$t(s) = iCtis)Ms) + Qt(s), (4.22) 
as 

where 

(4.23) 

Qt{s)^-^ E Y.^rn[ft^s{z)]as{al,{W{f))) {a.{d.Vz),W{g)} . 
zgAl zez 

Arguing as before, we arrive at the bound 

\\{atiWif)),Wig)}\\^ < \\{anWif)),Wig)}\\^ (4.24) 

+ E E r |Ini[/t-«(.)]|||K(a.Fz),lF(g)}|Uds. 

Inserting (|4.8p into (|4.24p and using the harmonic bounds from Corollary [10] with = /is, 
we find that 

||{a*(W(/)),lF(5)}||^<C||/||oo||5llooe"* ^ F^M^,y)) 

xex,yeY 

t I- _ (4.25) 

+ CII/II00 E EE^A.3«^,^)) / e*(*-^) / \d.V{r-Z)\\\{as{W{r-5z)),W{g)}\\^drds, 
Z(ZALzeZxeX ■'^ -^^^ 

with C as in (|4.12p . and we have set v = {fi^ + e)wh(M3 + ^) for notational convenience. 
After iterating (|4.25p m > 1 times, we find that 

m 

\\{at{W{f)),W{g)}\\^<C\\f\\^\\g\\^e'' ^ J] a„(x, y; t) + i?„,+i(t), (4.26) 

where 

ao{x,y;t)^F^,{d{x,y)), (4.27) 
SLi{x,y;t) = CtY, E / \\r-5z\\oo-\d^{r-,Z)\drF^,{d{x,zi))F^,{d[z2,y)) 

< Ct J2 F^^id{x,z^))F^,id{z2,y)) E / M-|W(r;Z)|dr 

zi,Z2<£Al ZcAl •'^^ 

Zi,Z2&Z 

< CCst J2 F,Ad{x,Zi)) F^,{d{zi,Z2)) F^,{d{z2,y)) 

Zi,Z2&Al 

< CC:,CltF^,{d{x,y)) , (4.28) 
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and in general, 

{cty 



an{x,y]t) 

III . 

X (d(a;,zi,i)) • i^s (^^(21,2, ^2,1)) • ' ' -^^3 {d{zn,2,v)) 

- E E "' E -P'm • -^A's ('^(^1.2,2:24)) •• •-F'm (C?(^n,2,y)) 

Z1,1:Z1,2GAi, Z2,1,Z2.2&!^L Z„,i,Z„.2eAi 

|r,|.|W(r,;Z,)|dr, 



< 



xR E 

J = l Z,<ZK-L 



E 



,1 ^ 

^1,1 7^1,2,^2, 1,^2,2,' ■■ ,ZTi,l,2Ti,2f=A 



Zi^i)) • f)_,3 (d(zi4,zi^2)) • ^)<3 {d{zi,2, Z2^i)) ■ ■ ■ F^^ (d{zn.2,v)) 



< -, Ft^3 (d{x,y)) , 



(4.29) 



for any 1 < n < m. As before, with t > fixed, the remainder term i?j„+i(<) converges to zero 
as m — > cx). Thus we have proven (|4.1ip as claimed. □ 



5 A Priori Solution Estimates 

In this section, we will prove a variety of a priori estimates which will be useful in our proofs of 
the main results. The underlying argument which facilitates most of the lemmas below is the 
well-known Gronwall inequality. We state and prove a version of this estimate which is tailored 
to the present work. A more general bound of this type appears, e.g. in [1]. 

Lemma 14 (Gronwall Inequality). Let u : R ^ C satisfy 

\u{t)\<a{t)+ I /(t,s) |u(s)|ds (5.1) 

J a 

for all t in [a,b]. If a is non-negative and non- decreasing and f is non-negative and continuous 
with /(■, s) nondecreasing for each fixed s G [a, b], then 



|u(OI <a(i)exp(^y f{t,s)dsj . (5.2) 

for all t in [a, b] . 

Proof. We prove (|5.2p pointwise. Let t^ £ [a, b] and observe that 

\u{t)\<a{to)^ I f{t^,s)\u{s)\ds, (5.3) 

J a 

holds for all t E [a, to]. Define 

m{t)^a{to)-^ f fito,s)\u{s)\ds. (5.4) 



Clearly, \u{t)\ < m{t) and the bound 

m\t)^ f{to,t)\u{t)\ < f{ta,t)m{t), (5.5) 
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readliy implies 

\u{t)\ < m{t) < m(0) exp i^J /(to, s) ds ] , 
for any t G [a, to]- Taking t — to, we have proven (|5.2p 



(5.6) 
□ 



Recall that our general, finite volume, multi-site Hamiltonian, H^^ : K, is of the form 



The applications we have in mind concern bounding the solutions of our Hamiltonian flows. 

imiltonia] 

E (5-7) 



ZcAz 



and we need a variety of assumptions on the perturbations Vz to prove our estimates. 

We begin with a basic proof of boundedness for the flow $4 : X\j^ corresponding to 

(|5.7p . As is demonstrated in [9], boundedness follows if the perturbation is dominated by the 
harmonic part. For the sake of completeness, we include this argument here. 

We assume the perturbation in (|5.7p above satisfles: There exist numbers Ci > 0, Ci > 0, 
and Hi > such that 



(5.8) 



for each a; G and any x G .Ya^, . 

Lemma 15. Fix L > 1 and let $f denote the flow corresponding to the Hamiltonian H^'- 
defined in \5. 7| ) above. If the perturbation satisfies \5.^) described above, then for any x G Aa^^ 
the components of the flow $t(x) = {(cix{t),Px(t))}xe!^L satisfy 



sup max(|gf;i(t)|, \px{t) \ ) < Ki exp{K2t), 



where 



and 



Ki ^Ki{x)^J sup (p2(o) + ^2(0) + Ci 
xeA 



j=l xGAl 



(5.9) 
(5.10) 

(5.11) 



Proof. Fix L > 1, take x G A^, and choose x G Xa^,- Consider the function deflned by setting 

Exit) = pl{t) + qlit) + Ci , (5.12) 
where Ci > is the number appearing in (j5.8p . From Hamilton's equations, we have that 

Exit) = 2pxit)pxit) + 2qxit)qxit) 



-4 L;2^2EAj - 1 Pxit)qxit) + '^Pxit) E W + (lx-e,{t)) 



-2pxit) E [dxVz] ($t(x)) , 



and therefore, 



Exit) < E ^-.y^yW = (AEit))^ 
yeAi, 



(5.13) 



(5.14) 
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where the |Ai| x |Ai| matrix A = {A^^y) is given by 
2' " 



A. 



x,y-) 2Xj + Ci Ff,,{l) ify = x±ej, 

Ci F^^ {d{x,y)) otherwise, 



(5.15) 



and {AE{t))x is the x-th component of this vector. Denote by E the vector- valued fimction 
whose components are and equip W with the sup- norm || • I loo- With ([5?T4ll it is clear that 



and therefore, 

mm^ < \\Eio)i 



m < \\AE{t)\\ 



E{s) ds < WEmic 



\AE{s) 



Letting u{t) — \\E{t)\\oo, Lemma [Ml implies 

max(|g.(i)|Mp.(i)P)<||i?(t)|l^ < ||i?(0)|U exp (P|U t) 
from which (|5.9p is clear. 



(5.16) 

(5.17) 

(5.18) 
□ 



As will become clear in the proof of Lemma [T7] below, the main quantities of interest for 
us are the derivatives of the components of the flow with respect to the initial conditions. The 
next lemma provides explicit estimates on these functions. To prove it we need the following 
additional assumption on our perturbation. Assume there exists constants C2 > and ^2 > 
for which: given any L > 1 and any pair a;, y G Aj^, the bound 



\[d.dyVz]i^)\ < C2F^Ad{x,y)) , 



(5.19) 



holds for all points x G Xf^j^ . 

Lemma 16. Fix L > 1 and let denote the flow corresponding to the Hamiltonian H^^ 
defined in |5. 7| ) above. If the perturbation satisfies I15.8\) and i5.19\) . then for any x G X\,^ the 
components of the flow $t(x) = {iqxit),Pxit))}xi^AL satisfy 



sup max 



dqx{t) 



dqy{0) 



dqx{t) 



dpy{Q) 



< max(l,2i) exp^Kt"^), 



(5.20) 



and 



sup max 

x,y£Az, 



dpx{t) 



dqyiO) 



dpx{t) 



dpyiO) 



<l + tjx + 2^Ajj max(l,2t) exp(Xt2), (5. 21) 



where 



K = 2^2 



i=l x^Kl 



(5.22) 
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Proof. We begin with a proof of (|5.20p . Using Hamilton's equations, wc find that 



Qxit) - fe(0) - 2tp^{0) 



Px{r)dr ds 



(5.23) 



f r-t 

-2 ^ / (i-s)[a,Fz]($.(x)) ds 

. "'0 



and therefore 

dqyiO) 



S4y)-A I + 2^EA, j 
+ 4^ A, f\t-s) 



dqy{0) 
dqx+ej (s) dq^^ej (s) ' 



(5.24) 



zGAl ZcA 

More succinctly, we have found that 



dqyiO) dqyiO) 

,c A , A , "'0 



ds 

'As) 



,(0) 



ds. 



dqxit) 



.(0) 



<4(y)+ E 

.^A "'0 



zeAi, 



dqz{s) 



where the jA^j x jA^j matrix A = {A^ z) has entries 



A^ 



dqyiO) 



4U^' + 2E;=iA, +2C2F^,(0) if. 



ds 



(5.25) 



4A, + 2C2F^,(1) 
2C2F^, id{x, z)) 



If, for each fixed y G A^, we denote by dyq 



(5.26) 

^l^^l the vector- valued function whose 



if z = X ± Cj , 
otherwise. 



dqyiO) 



, then (|5.25p implies 



components are given by 



Again letting u{t) — \\dyq{t)\\oo, Lemma [14] yields the estimate 



Rq(t)|U<l + / {t-s)\\Adyq{s)\\^ds. 



||5j,g(t)|L < exp 



Quite similarly, the bound 

'qx{t) 



dpyiO) 



< 2t5M+ E / (^-«)^ 



dqy{Q) 



ds . 



(5.27) 



(5.28) 



(5.29) 



follows from (|5.23[) with the same matrix A. This proves (|5.20p . 
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The bound for Px{t) follows from (|5.20p . In fact, it is easy to see that 

Px{t) = px{0)-2(u;^ + 2j2x}j q4s)ds (5.30) 



+ 2^A, / (fe+e,(s) + fe-e,(s)) ds - / [dxVz]{^s{x))ds. 

Using dSSni), (jOT]) readily follows. □ 

Lemma 17. Let X,Y C Z'^ be finite sets and take Lq > 1 large enough so that X,Y C A^^ . For 
any L > Lq and t G M, denote by the dynamics corresponding to ( [5. 7\ l. If the perturbation 
satisfies i5.8\) and i5.19\) . then there exist positive numbers Ki and K2, both independent of L, 
for which: given any functions f : X C and g :Y —> C, the bound 

\\{ai{W{f)),W{g)]\\^ < ifi|X||r|||/|U||.g||oo exp(i^2i'), (5.31) 
holds for all t CzR. 

Proof. We first fix L > Lq as in the statement of the lemma and prove the estimate on A^,. 
In this case, we suppress the dependence of most quantities on L to ease notation. Now, recall 
that for any fixed point x, 

[{at{Wif)),Wig)}] (x) = ^ A [a,{Wif))] (x) • A [14/(5)] (x) 

-E^Mm/))](x).^[I4^(.)](x). 

Since 

A [W{g)] (x) = zim [5(2;)] [W{9)] (x) and A [w{g)] (x) = iRe [g{v)] [W{g)] (x) , (5.33) 

OPy Oqy 

the sums in ()5.32p above are only over those y in the support of g. The derivative of the 
time-evolved quantities may also be calculated, e.g., 

A [a,{W{f))] (x) = * [ ^ Re[/(x)] ■ + Im[/(x)] ■ .[a,{W {f))] (x) • (5.34) 

oqy oqy oqy J 

Clearly, the sum in (|5.34p is only over those x in the support of /, and a similar formula holds 
for the derivative with respect to py . Thus, 

irr i^xrifw Mi^ M ^ AWf ( ^9^(0 dp^t) dq^t) dp^t) \ 

\[{at{W{f)),W{g)}]{y.)\ < 4 H^IU 2^ max — , — , — , — 

xex,y<£Y ^ '^^ '^y ' 

Using Lemma fTHl (|5.3ip immediately follows. □ 
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